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Abstract: In this paper we study the Sudakov form factor in = 4 super Yang-Mills 
theory to the three-loop order. The latter is expressed in terms of planar and non-planar 
loop integrals. We show that it is possible to choose a representation in which each loop 
integral has uniform transcendentality. We verify analytically the expected exponentiation 
of the infrared divergences with the correct values of the three-loop cusp and collinear 
anomalous dimensions in dimensional regularisation. We find that the form factor in A/" = 4 
super Yang-Mills can be related to the leading transcendentality part of the quark and 
gluon form factors in QCD. We also study the ultraviolet properties of the form factor in 
D > A dimensions, and find unexpected cancellations, resulting in an improved ultraviolet 
behaviour. 



Keywords: 



Bupersymmetric gauge theory, NLO Computations . 



Contents 



|l]. Introduction ^ 

|2|. Form factor to two loops ^ 

|3|. Momentum routing invariances of integrals 

^ Form factor to three loops from unitarity cuts [To| 

[4.1| One-loop form factor from unitarity cuts 11 



4.2 Two-loop form factor from unitarity cuts 12 



4.3| Three-loop form factor from unitarity cuts 13 

|5|. Final result for the form factor at three loops 15 

|6|. Logarithm of the form factor 17 

|7|. Ultraviolet divergences in higher dimensions 18 

|8|. Discussion and conclusion 21 

|A] . Explicit results of integrals 23 

1^. Form factor in terms of master integrals 28 

10. Four-point amplitude to two loops 29 



1. Introduction 

In this paper we study the Sudakov form factor in = 4 super Yang-Mills (SYM) with 
gauge group SU{N). Following van Neerven [1], we study the vacuum expectation value of 
an operator built from two scalars, inserted into two on-shell states. The operator belongs 
to the stress-energy supermultiplet, which contains the conserved currents of = 4 SYM, 
and has zero anomalous dimension. Together with the vanishing (3 function of = 4 SYM 
this means that the form factor is ultraviolet (UV) finite in four dimensions. Therefore 
only infrared (IR) divergences associated to the on-shell states appear, which we regularise 
using dimensional regularisation. 

Generalisations of the Sudakov form factor to the case of different composite operators, 
and more external on-shell legs, have been discussed recently in refs. [2-5]. Form factors 
have also been studied within the AdS/CFT correspondence in the dual AdS description, 
see refs. [6,7]. Here we will focus on the perturbative expansion of the form factor of ref. [1]. 
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Form factors are closely related to scattering amplitudes. For example, planar ampli- 
tudes can be factorised into an infrared divergent part, given by a product of form factors, 
and an infrared finite remainder (a 'hard' function in QCD terminology), see e.g. ref. [8] and 
references therein. The infrared divergent part exponentiates and has a simple universal 
form. In fact, for four- and five-point scattering amplitudes the exponentiation property 
of the divergent part carries over to the finite part as well [8,9]. This is a consequence of 
a hidden dual conformal symmetry of planar scattering amplitudes. The latter relates the 
finite part to the infrared divergent part through a Ward identity [10, 11]. The relation to 
form factors makes it possible to give an operator definition of the finite remainder. The 
scheme independence of the latter was recently checked in a two-loop computation using 
dimensional and massive regularisations [12]. 

Scattering amplitudes in = 4 SYM have many special properties, and it is interesting 
to ask how much of this simplicity carries over to the form factors. For both the planar four- 
particle amplitude and the form factor, the general form of the result is known in principle. 
For the former, this is due to dual conformal symmetry, and for the latter it is due to the 
exponentiation of infrared divergences. However it is quite non-trivial to obtain these a 
priori known results from explicit perturbative calculations, evaluating loop integrals. The 
simplicity of the final results suggests that there should be more structure hidden in the 
loop integral expressions, and by studying them further one might gain insights into better 
ways of evaluating them, which is of more general interest. 

One might expect that the evaluation of form factors should be simpler than that of 
scattering amplitudes, as the former have a trivial scale dependence only, whereas the latter 
are functions of ratios of Mandelstam variables, e.g. s/t in the four-point case. Given this, 
it is somewhat surprising that less is known about the loop expansion of form factors in 
AT = 4 SYM than about scattering amplitudes. For example, while the planar four-point 
amplitude was evaluated to the four-loop order (in part numerically) [13-15], the form 
factor has only been computed to the two- loop order in ref. [1] , in a calculation that dates 
back to 1986. In the present paper, we extend the calculation of ref. [1] to three loops, 
and study which of the properties that have been observed for scattering amplitudes are 
present. 

One fact which makes form factors technically challenging compared to planar ampli- 
tudes, however also more interesting, is the following. At leading order in the 't Hooft 
limit N ^ oo, where the coupling A = g^N is kept fixed, both planar as well as non-planar 
integrals appear in the form factor. This is easily understood by the fact that the operator 
insertion is a colour-singlet. It is interesting to note that the non-planar diagrams appear- 
ing in the form factor are related, through the unitarity technique, to a priori subleading 
double trace terms in the four-particle scattering amplitude. Therefore, the form factor at 
leading order in N contains information about non-planar corrections to the four-particle 
amplitude. The first non-planar diagram, the crossed ladder, appears at the two-loop level. 
At three loops, we find five different non-planar diagrams that contribute, i.e. that have 
non- vanishing coefficient. 

It is an observed, albeit unproven fact that results for scattering amplitudes in A^ = 4 
super Yang- Mills have uniform transcendentality (UT), i.e. can be expressed as linear 
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combinations of polylogarithmic functions of uniform degree 2L, where L is the loop order, 
with constant coefficients. In e-expansions of dimensionally regularised quantities which 
depend only on a single scale, the coefficients of the Laurent expansion in e are real constants 
which are in general of increasing transcendentality in the Riemann ^-function. In this 
context uniform transcendentality refers to homogeneity in the degree of transcendentality 
(DT), where the latter is defined as 



In the planar case, the property of UT is even true for individual loop integrals, at least 
when they are expressed in an appropriate basis of dual conformal integrals [16,17]. Inci- 
dentally, this also has practical advantages, as these integrals are easier to evaluate [17, 18] 
than those in other representations. Dual conformal symmetry is only expected in the 
planar case, but what can be said about the transcendentality properties of non-planar 
integrals? At four points, the non-planar double ladder integral is not of uniform tran- 
scendentality. However, if defined with an appropriate loop-dependent numerator factor, 
it does have this property [19,20]. Changing to a basis involving the latter integral allows 
one to understand the UT property of four-point non-planar A/" = 4 SYM amplitudes [21] 
and J\f = 8 supergravity amplitudes [22,23]. It also raises the interesting question whether 
this is a generic feature. 

All planar and non-planar master integrals for form factors in dimensional rcgularisa- 
tion at three loops are known from the computation of the form factor in QCD [24-33] , and 
some of them have UT, while others do not. It has been observed that some of the inte- 
grals do have UT if they are defined with certain (loop-dependent) numerator factors [19]. 
The latter resemble the numerator factors required by dual conformal symmetry in the 
planar case [34]. In this paper, we find similar numerator factors for all topologies with 
7, 8, 9 propagators, such that the integrals have UT. Moreover, we find that the complete 
three-loop form factor can be written solely in terms of UT integrals. 

Finding a representation that has this property required using certain identities for 
non-planar form factor integrals that are based on reparametrisation invar iances, which we 
found as a by-product of our analysis. They generalise an identity found by Davydychev 
and Usyukina [35]. 

As was already mentioned, in = 4 SYM, scattering amplitudes and the form factors 
studied here are UV finite in four dimensions. It is interesting to ask in what dimension, 
called critical dimension Dc, they first develop UV divergences. This question is of theoret- 
ical interest in the context of the discussion of possible finiteness of A/" = 8 supergravity, see 
e.g. [36-38] and references therein. More practically, bounds on the critical UV dimension 
at a given loop order can also be a useful cross-check of computations, or constrain the 
types of loop integrals that can appear. Ultraviolet power counting, based on the existence 
of M = 3 off-shell superspace [39], provides a lower bound for the critical dimension. We 
analyse the UV properties of the form factor to three loops and find that at each loop 



DT{r) 
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DT{x) + DT{y) . 
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order, the critical dimension is Dc = 6. This is consistent with the bound obtained from 
superspace power counting. We find that the latter bound is saturated at two loops, while 
it is too conservative at three loops, where the ultraviolet behaviour is better than sug- 
gested by the bound. This is the result of a cancellation between different loop integrals. 
We find a representation where the UV behaviour is manifest. 

This paper is organised as follows. We review the known expression for the form factor 
to two loops in section ^. We then discuss identities for non-planar integrals to three loops 
in section ^. In section ^, using the unitarity-based method, we derive an expression for 
the three- loop form factor in terms of loop integrals. We then evaluate the latter in section 
|5| and verify the exponentiation of infrared divergences in section ^. We then analyse the 
ultraviolet properties of the form factor to three loops in section 0. We conclude in section 
|8[ There are several appendices. Appendix |A| contains the analytic expressions of the e 
expansion of the integrals used in the paper, while appendix]^ contains the expression of the 
form factor in terms of conventionally used master integrals. Finally, appendix reviews 
the on-shell four-point amplitude to two loops that is used in the unitarity calculation in 
the main text. 

2. Form factor to two loops 

In order to define the scalar form factor in A/" = 4 SYM, we start by introducing the bilinear 
operator 

O = Tr((/<i2(/.i2) , (2.1) 

where the scalars (pAB are in the representation 6 of 5?7(4), and (jy^s = ^'ab'^o,, with Ta 
being the generators of SU{N) in the fundamental representation, normalised according 
to Tr(r'^r'') = S'^K This operator is a particular component of the stress-energy super- 
multiplet of A/" = 4 SYM, and has zero anomalous dimension. We then define the form 
factor as the vacuum expectation value of O inserted into two on-shell states in the adjoint 
representation, 

Ts = mM44{P2)0), (2.2) 

with the convention that momentum is outgoing. 

Since O is a colour singlet, the form factor must be proportional to Tr (T"T^), 

= TY(T"r^)F5 . (2.3) 

We work in dimensional regularisation with L) = 4 — 2e dimensions in order to regulate IR 
divergences associated with the on-shell legs. We write the form factor as an expansion in 
the 't Hooft coupling [8] 



87r2 

according to 



{Any e-'^^ , (2.4) 



Fs = l + ax' F^^^ + x^' F^ + x^' F^ + 0{a^) . (2.5) 
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We normalized the tree-level contribution to unity and introduced 

X = f . , (2.6) 
—q'^ — IT] 

with the infinitesimal quantity 77 > 0. 

We remark that the dependence on the number of colours N in equation (2.5) is exact. 
In order to see this, let us show that the three-loop contribution to the form factor must 
be proportional to (a similar analysis trivially holds at one and two loops). 

The reasoning is as follows. Imagine a generic Feynman diagram contributing to J-s- 
Without loss of generality, suppose that it is built from three-point vertices, whose colour 
dependence is given by the structure constants /"i^a"^. For each internal line, there is a 
sum over adjoint colour indices, with the result being proportional to Tr(T"T^), as stated 
in equation ( ^.31) . Our goal is to determine the proportionality factor. In order to do this, 
it is convenient to sum also over the free indices a and 6, 

^ j«6Tr(r''r*) = - 1 . (2.7) 

a,b 

We can then represent each Feynman diagram as a circle with inscribed lines. There are 
three inequivalent structures that can appear, 

jabg jbcg jcdh jedi J^f^ jfo-^ 

— jabg j^bch jcdg jdei jeif jfha ^ C^.S) 
^ jabg jbch jcdi jdeg j^^f jf^'^ 

which correspond to the case of zero, one, or two intersections of the inscribed lines, 
respectively. Sums over repeated indices are implicit. In order to carry out the sums, it is 
convenient to write the structure constants as 



= -i/V2 (TV[T"r^r'=] - Tr [T^T^T^]] . (2.9) 
Using the SU{N) Fierz identities, 

Y^Tt:{AT'')Tt:{BT'') = Tt{AB) - 1/iV Tr(yl) Tr(fi) , (2.10) 

a 

^Tt{AT''BT°') = Tt{A)Tt:{B) - 1/NTiiAB), (2.11) 



one easily finds 



(iV2-l)iV3, B = -^{N'^ - 1)N^ , C = 0. (2.12) 



Taking into account equation ( |2.7D , we see that Fs at three loops is proportional to A^^, 
as claimed. 

Note that beginning from four loops there can be more than one colour structure, and 
in particular the quartic Casimir can appear. An explicit example of this is the four-loop 
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Di El E2 

Figure 1: Diagrams tliat contribute to ttie one- loop and two- loop form factors J^^^' and J-'g^ in 
M = A SYM. All internal lines are massless. The incoming momentum is q = Pi + P2, outgoing 
lines are massless and on-shell, i. e. Pi — P2 — 0- All diagrams displayed have unit numerator and 
exhibit uniform transcendentality (UT) in their Laurent expansion in e = (4 — D)/2. 



contribution to the QCD /3 function [40]. An interesting related question has to do with 
the colour dependence of infrared divergences in gauge theories, see e.g. [41], and references 
therein. 

The form factor to two loops was computed a long time ago [1]. It contains as building 
blocks the diagrams displayed in Fig. || and reads 



= l + ax'R,-2 P + a2 x^' ■ [4 p + E^""^] + ©(a^) , (2.13) 

Re = 7 T • (2.14) 

2r(l-e) ^ ^ 

The expressions for Di, D^^^, Ei, and E^^^ are given explicitly in appendix ^ and result 



m 



/3IC5 _ 77r^\ ^4 / 9497r6 _ 49C| \ 
V 5 36 y ^ V 120960 18 J 

% _ 3l7r^C5 127CA ^6 /^ 497r^C| _ 217C3C5 185937r8 \ 



'e2 ^ 12 ^ 3 1440 ' ' I + e 



5/ 3297r ^ ^ 

+^ V 4320 60~ + 7 J +^ 216 15 + 9676800 ) 

+0{e') , (2.15) 



fP = i?2.[4^fP + £;fP] 

^ I ^237r2C3 TlCsN 2 ^901C| 2577r6 \ 

2e^ 24e2 12e 240 A 72 20 7 V 36 6720 J 

3 / 129l7r4C3 _ 3137r^C5 3169C7 \ 

V 1440 120 14 J 
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/ 845C3C5 15477r2c? 504197r«\ 

-66C5 3H ^ — + ]+0(e^). (2.16) 

V 6 216 518400 / 

The multiple zeta values Cmi,...,mk defined by (see e.g. [42] and references therein) 

The numerical values of the transcendental constants up to weight eight are: 

Ca = 1.2020569031595942854 . . . , Cs = 1.0369277551433699263 . . . , 
Ct = 1.0083492773819228268 . . . , C5,3 = 0.037707672984847544011 .... 

3. Momentum routing invariances of integrals 

Before we proceed to calculate the A/" = 4 SYM form factor to three loops via unitarity 
cuts, we want to investigate some of the occurring topologies more closely. In particular, 
we will derive identities that relate integrals without uniform transcendentality (UT) to 
integrals that do have this property. Since the diagrams that we will obtain from the 
unitarity method do not individually have UT, the following relations will be very useful 
later on for switching to an integral basis for the form factor in which each building block 
has UT. 

We start with topology F^, see Fig. |^. We label its incoming momentum with q = 
Pi + P2, and the outgoing ones with pi and p2, respectively. The latter are massless and 
on-shell, i.e. pf = P2 = 0. The topology can be parametrised according to 

{ki - k2 , ki - ks , ki ~ k2 - k-s , k2 , h , ki - q , k2 - q , h - q , k2 - pi} , (3.1) 

where ki are the loop momenta. It can be seen from Fig. § how the momenta are distributed 
among the lines of the diagram F^. It turns out that the following reparametrization of 
loop momenta, 

ki q + k2- ki 

k2 k2 

h ^ q-h , 

does not only leave the value of the integral invariant, but even its integrand. We can 
now apply this transformation to the integral F^ which carries the factor (^2 — ^3)^ as an 
irreducible scalar product in its numerator. This yields 

ik2 - ^ {k2 + k3- qf 

= kl + kl + {k2 - qf + (fca - qf - {k2 - ksf - q^ . (3.2) 

We can now solve this equation for {k2 — k^f and get the following relation between 
integrals, 

F3 = -^q'F^+Fai+Fs, (3.3) 
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Figure 2: Diagrammatic representation of Eq. (3^). The internal lines of all diagrams are massless. 
The incoming momentum is q ~ pi + P2, outgoing lines are massless and on-shell, i.e. Pi=P2 — 0- 
Diagrams with labels pa and pb on arrow lines have an irreducible scalar product {pa +Pb)^ in their 
numerator (diagrams that lack these labels have unit numerator) . The numbers in indicate the 
position of the entries in Eq. (3.1). Diagrams F3 and Fg have UT, contrary to F3* and Fai- 



which is diagrammatically shown in Fig. |2|. We have now decomposed the integral -F3, 
which does not have UT in its Laurent expansion, into two integrals (F3 and Fg) which 
indeed do have this property, and the auxiliary integral Fai, which again does not have 
homogeneous transcendental weight, but which will be cancelled later on. 

We can apply analogous steps to topology F|, see Fig. |3[ The topology can be 
parametrised according to 

{h , k2 , h , ki - k2 , h - ks , ki - q , ki - k2 - P2 , h - q , k2 - pi} , (3.4) 

and the distribution of the momenta among the lines can be seen from Fig. |3|. The integrand 
remains invariant under 

ki ^ q — ki 
h ^ Pi - h 
h ^ q-h , 

We now apply this transformation to the numerator (ki — pi)^ of the integral F4. This 
yields 

{h-pif {h-P2? 

= kl + {ki-qf -{ki-pif -q^ . (3.5) 
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Fa2 



Figure 3: Diagrammatic representation of Eq. (3.6). All symbols have the same meaning as in 
Fig. ^. The numbers in indicate the position of the entries in Eq. (3.4). Diagram F4 has UT, 
contrary to Fl and Fa2- 



We can now solve this equation for {ki — pi)^ and get 

F^ = -^q^F*+Fa2, (3.6) 

which is diagrammatically shown in Fig. ^. Again we decomposed the non-homogeneous 
integral F^ into the homogeneous integral F4 and yet another non-homogeneous auxiliary 
integral {Fa2) which will be cancelled later on. 

We can also decompose the topology F^, see Fig. ^ In this case we cannot find a 
relation between integrals which is based on a momentum routing invariance, but a relation 
which is simply based on momentum conservation. The topology can be parametrised 
according to 

{ki - k2 , ki - k3 , ki - k2 - k3 , k2 , ks , ki - q , k2 - q , h - pi , ks - pi} , (3.7) 

and we refer to Fig. |^ for their distributions among the lines. From momentum conservation 
we get 

{k2 - pif = {ki - k2f -kl + kl + {ki - pif - {ki -k2- pif , (3.8) 
which results in 

F* = Fal + Fa2 +F^-F^-Fq. (3.9) 

Hence we decomposed F^ into the homogeneous- weight diagrams -F5, Fq, and Fg, as well 
as the same non-homogeneous diagrams F^i, and Fa2 which already appeared above. 

We see from Eqs. (p.3|), (|3.6D , and (|3.9D that only two auxiliary diagrams of non- 
homogeneous weight, namely Fai, and Fa2 appear in all these relations. It turns out that 
the coefficients obtained from unitarity are precisely such that these integrals cancel in the 
expression for the form factor. 

We checked all relations between integrals also at the level of their integration-by- 
parts (IBP) reduction [43] to master integrals using the implementation of the Laporta 
algorithm [44] in the REDUZE [45] code. We find that all relations obtained from mo- 
mentum routing invariance in this section can actually be reproduced from solving IBP 
relations, which is a priori not guaranteed for a general Feynman integral topology. The 
e-expansions of all integrals can be found in appendix |A[ 
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Figure 4: Diagrammatic representation of Eq. (3.9). All symbols have the same meaning as in 
Fig. ||. The nmxibers in indicate the position of the entries in Eq. ([3.7]). Diagrams F5, Fg, and 
Fg have UT, contrary to F5*, Fai, and Fa2- 



4. Form factor to three loops from unitarity cuts 



Here we use unitarity cuts to derive an expression for the three-loop form factor in terms 
of the integrals discussed in the previous section. We will compute the form factor in a 
perturbative expansion in the Yang-Mills coupling g, and denote the contribution at order 
, , , by J^g^^,F^ loop ^ j~3- loop ^ respectively, and similarly for Fs- Note 

that this notation, convenient for the unitarity calculations, differs from the one used in 
Eq. (P). 

The essential features of the unitarity-based method [46,47] that we are going to use are 
reviewed in the recent paper [48]. We will employ two-particle cuts, as well as generalised 
cuts. The two-particle cuts are very easy to evaluate, and we show an explicit example 
below. 

In order to evaluate more complicated cuts, with many intermediate state sums to 
be carried out, it is extremely useful to employ a formalism that makes supersymmetry 
manifest. This can be done by arranging the on-shell states of = 4 SYM into an on-shell 
supermultiplet [49]. The main advantage is that intermediate state sums appearing in the 
cuts become simple Grassmann integrals that can be carried out trivially [50-52]. In this 
way, it is easy to obain compact analytical expressions for the cuts. 

We follow the notations for unitarity cuts of ref . [53] . We start by reviewing the one- 
and two-loop cases as examples. 
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2 — loop 





-pl-loop A\~^^"^ 

(36) 




^2 — I oop 



(3c) 








Figure 5: Two-particle cuts of form factors up to tlrrce loops. 



4.1 One-loop form factor from unitarity cuts 

As a simple warmup exercise, we rederive the one-loop result from unitarity cuts, see also 
ref. [2]. Let us compute the two-particle cut (la) shown in Fig. |5[ It is given by 



I— loop 



S 



cut(la) 



E 

Pl,P2 



^ -ntreer n 



(4.1) 



where £i and £2 are the momenta of the cut legs, and the sum runs over all possible particles 
across the cut. We may use the on-shell condition £f = £| = in the integrand (but not on 
the cut propagators), since any terms proportional to such numerator factors would vanish 
in the cut. The four-particle tree ampliutde .4.4''^'^ (£2, ^15^1,^2) is given in appendix]^. We 
use the convention that all momenta are defined as outgoing. 

When computing the cut of a form factor (as opposed to a colour-ordered amplitude), 
one has to be careful about the overall normalisation, since the possible exchange of external 
legs pi and p2 leads to a factor of 2 in the cuts. When comparing cuts of the form factor 
to cuts of integrals, this factor cancels out. In the following we count such contributions 
only once. 

The two-particle cuts are particularly simple to evaluate. With our choice of external 
states, only scalars can appear as intermediate particles, and we therefore do not need to 
use the spinor helicity formalism. The tree-level form factor is simply given by 



j"|r"<^(-4,-^2) = Tr(r"r^ 



(4.2) 



The necessary four-particle amplitudes are given in appendix j^. The colour algebra across 
the cut is carried out using the SU(N) Fierz identities, see eqs. ( p. 10 ) and ( p. 11 ). It is 
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Figure 6: Diagrams that do not contribute to the form factor at two (£^3) and three loops (F7 
and i^io)j respectively. They have worse UV properties compared to the integrals that do appear 
in the form factor. The labels pa and pb on £7 indicate an irreducible scalar product {pa + Pb)'^ in 
its numerator. The other two diagrams have unit numerator. 



easy to see that (4.1) becomes 



J c 



c/V^'^^si2Tr(T"T'' 



cut (la) 



d^k i i 

d^'k 



-Ig^fi^^NsuTriT'^T") [ 

J «(27r 



(Pi+h) 
1 



+ 



{P2 + 



(2^)^P(fc+pi)2(fc-p2)2 



cut (la) 



ej=q=o ' 
(4.3) 



where Sij := {pi + Pj)"^ , and where we have identified the cut of the one-loop form factor 
with the cut of the one-loop triangle integral Di, see Fig. ||, 

d^k 1 



Di 



(4.4) 



i(2^)^ k'^{k+pif{k-p2Y ' 

We can now argue that this result is exact, i.e. that we can remove the "cut (la)" in 
Eq. ( [4.31 ). In order to do that, we have to make sure that no terms with vanishing cuts are 
missed. Such terms having no cuts in four dimensions can be detected in D dimensions. 
The two-particle cut calculation we just presented would have gone through unchanged in 
D dimensions, since all required amplitudes were those of scalars, and no spinor helicity 
identities intrinsic to four dimensions were used. A similar argument was given in ref. [53]. 
Therefore we conclude that in D dimensions, 



rpl — loop 



g'Nfi''{-q')2Di . 

4.2 Two-loop form factor from unitarity cuts 

We recall that at two loops, the result for the form factor is given by [1], 



(4.5) 



(4.6) 



where the planar and non-planar ladder diagrams Ei and E2 are shown in Fig. ||. 

Let us now understand this result from unitarity cuts. The unitarity cut (2b) of Fig. ^ 
detects the presence of the planar integral Ei only. The calculation is identical to that 
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of the one-loop case, with the exception that the one-loop form factor as opposed to the 
tree-level form factor is inserted on the l.h.s. of the cut. 

The unitarity cut (2a) of Fig. ^ reveals a new feature, that was already mentioned 
in the introduction. On the r.h.s. of the cut we now insert the full one-loop four-point 
amplitude ^4"'°°^, given explicitly in Eq. (C.2), which in addition to single trace terms also 



contains double trace terms. The latter would ordinarily be subleading in the expansion of 
powers of N, e.g. when computing a four-point amplitude at leading colour using unitarity 
cuts. Here, however the colour algebra gives rise to another factor of for those terms, 
so that they can contribute to the form factor at the same order as the single trace terms. 
This explains why the non-planar integral E2 can appear in the form factor. 

In principle, new terms could appear in the three-particle cut, but this is not the case. 
For example, the diagram shown in Fig. ^ has no two-particle cuts. The absence of this 
diagram can be understood by the fact that it has worse UV properties compared to Ei 
and E2, as we discuss in section ^. For the same reason, diagrams F7 and Fiq from Fig. |6|, 
the latter of which at has no three-particle cuts, will not contribute to the form factor at 
three loops, as we will see below. 

We have also evaluated the three-particle and a generalised cut, with the result being in 
perfect agreement with Eq. (|4.6| ) . We found it useful to employ a manifestly supersymmetric 
version of the unitarity method [50]. The necessary tree-level amplitudes for the local 
operator of Eq. ( p.l| ) inserted into three on-shell states were computed in refs. [2,3]. The 
analytical calculation is straightforward to perform. We refrain from presenting the details 
since it would require introducing spinor helicity and superspace. We refer the interested 
reader to refs. [48,50] for related instructive examples. 

4.3 Three-loop form factor from unitarity cuts 

We again begin by studying two-particle cuts, which are shown in the second line of Fig. |5|. 
Again, all results for the form factors and four-point amplitudes appearing in the unitarity 
cuts are explicitly known, with the result for the four-point amplitudes summarized in 
appendix 

When evaluating the cuts, one has a certain freedom in rewriting the answer to a given 
cut due to the on-shell conditions. Of course, eventually such ambiguities are fixed by the 
requirement that the answer must satisfy all cuts. In order to find such an expression that 
manifestly satisfies all cuts it is very useful to have an idea about the kind of integrals that 
should appear in the answer. We expect that the form factor can be expressed in terms of 
the integrals that have UT that were discussed in section |3[ This turns out to be a very 
useful guiding principle. 

The calculation is completely analogous to that at one and two loops. Let us start 
with the simplest cut (3c) from Fig. ^. It is given by 



-j~3—loop 



cut(3c) 



The evaluation of the cut is exactly as that considered at one loop, with the difference that 
we now insert the two-loop expression for the form factor into the cut, as opposed to the 



- 13 - 



tree-level one. One immediately finds 



3— loop 



cut(3c) 



^6^6. N^-qy[8Fi+2F*] 



cut(3c) 



(4.8) 



where Fi is the three-loop ladder integral shown in Fig. ^, and F^ is related to -F3 in the 
same figure via the identity (3.3). In fact, we know from section ^ that F| does not have 
uniform transcendentality. Since we do expect the final result to have this property, use 
Eq. (|3.3|) to eliminate F3. When doing so, we note that the contribution of Fai in that 
equation drops out on the cut (3c), and we have 



jpZ—loop 



cut (3c) 



cut(3c) 



(4.9) 



i.e. we have succeeded in writing the two-particle cut (3c) in terms of integrals having UT 
only. 

Similarly, one can show that the two-particle cut (3b) of Fig. ^ can be written as 



Tp3—loop 



cut(3b) 



/ ^3(_^2)2^g(_^2)^^ ^4^^j 



cut(3b) 



(4.10) 



This confirms the coefficient of Fi, and introduces a new integral F4, invisible to cut (3c). 

Finally, the most interesting two-particle cut is (3a), as it uses the double trace terms 
present in see appendix Q Using the identities derived in section |3|, we find 



^3— loop 



cut(3a) 



IJ.^' {-q^f[8{-q'^) Fi-2F2+AFs+AFi-AF5-AFe\ . (4.11 



cut (3a) 



Comparing equations ( [4.9| ),(4.10), and ( [4.11 ) with each other, we see that they are 
manifestly consistent with each other, which suggests that we are indeed working with an 
appropriate integral basis to describe this problem. We find that the following expression 
is in agreement with all two-particle cuts. 



rp3 — loop 



2— part, cut 



gG ^6. ^3(_^2)2^8(_^2)^^_2^2+4F3+4F4-4F5-4F6-4F8 



2— part, cut 

(4.12) 

It is quite remarkable that to three loops the coefficients of all integrals are small integer 
numbers. 

We could proceed by evaluating three- and four-particle cuts, but we find it technically 
simpler to study generalised cuts. To begin with, we perform a cross-check on the two- 
particle cut calculation above by evaluating maximal cuts where nine propagators are cut. 
We find perfect agreement between the two calculations. Next, we release one cut constraint 
to detect integrals having only eight propagators. There are several ways in which this can 
be done. For example, cutting all eight propagators present in integral Fg detects this 
integral, as well as integrals F5 and Fq. Another eight-propagator cut detects integrals 
F2,F5,Fq and Fj. The latter integral (see Fig. ^ turns out to have coefficient zero, i.e. it 
does not appear. 

We again find perfect agreement with the contributions already known from the two- 
particle cuts, and find further contributions not having any two-particle cuts, like Fg. The 
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following expression satisfies all cuts that we have evaluated, 



(4.13) 



We will now argue that Eq. (4.13) is the complete result for the three-loop form factor. In 
fact, potential corrections to equation ( f4.13| ) can come only from seven-propagator integrals 
that have vanishing two-particle cuts. An example of such an integral is Fio shown in Fig. |6|. 
As we will see in section |^, the appearance of such integrals is highly unlikely due to their 
bad UV behaviour, violating a bound based on supersymmetry power counting. 



Moreover, in section |6|, we will perform an even more stringent check on Eq. ( 4.13| ) by 



verifying the correct exponentiation of infrared divergences. In particular, this means that 



any potentially missing terms in equation ( 4.13 ) would have to be IR and UV finite, and 



vanish in all unitarity cuts that we considered. 

5. Final result for the form factor at three loops 

In the previous section we obtained the extension of Eq. ( |2.13| ) to three loops, 

Fs = 1 + g2 AT . (_^2) .2Di+g^N'' • {-q^f • [4 Fi + F2] 

+ / A^=^ ■ {-q^f ■ [8 (-g2) Fi - 2 F2 + 4 F3 + 4 F4 - 4 F5 - 4 Fg - 4 Fg + 2 Fg] 
+ 0(/) 

= l + ax'R,-2 + x"' RI ■ [4 F^^ + Ff P] 

W x^' Rl ■ [8 Ff P - 2 Ff P + 4 F|"P + 4 Ff p - 4 Ff p - 4 Ff p - 4 F|"p + 2 Ff P] 
+Oia^) . (5.1) 



The expressions for Fj, and F^^^ are again given in appendix^. All diagrams are displayed 
in Fig. 0. This yields 

F^^) = i?? • [8 Ff P - 2 Ff P + 4 F|"P + 4 Ff p - 4 Ff p - 4 Ff p - 4 F|"P + 2 Ff P] 

1 IIC3 2477r^ 1 / 857r2C3 439(5 \ 
~ ~6e6 ^ 12^3 + 25920e2 7 V 432 60~ J 

883C3 2252371^ ^ / 478037r^C3 244971^(5 385579(7 



36 466560 V 51840 432 1008 

2/1549 22499C3C5 , 4967r^C| 118375998l7rS \ , 3^ 
^^^5,3 + 7838208000 ) + ^ " ^^"^^ 

We can make a very interesting observation here. For anomalous dimensions of twist 
two operators, there is a heuristic leading transcendentality principle [54-56], which relates 
the A/" = 4 SYM result to the leading transcendental part of the QCD result. We can 
investigate whether a similar property holds for the form factor. 

For the comparison, we specify the QCD quark and gluon form factor to a super- 
symmetric Yang-Mills theory containing a bosonic and fermionic degree of freedom in the 
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Figure 7: Diagrams of which the three-loop form factor Fg in Af = 4 SYM is built. All internal 
lines are masslcss. The incoming momentum is q = pi +P2, outgoing lines are massless and on-shell, 
i.e. Pi = P2 = 0. Diagrams with labels pa and pt, on arrow lines have an irreducible scalar product 
(j>a + Pb)'^ in their numerator (diagrams that lack these labels have unit numerator). All diagrams 
displayed exhibit uniform transcendentality (UT) in their Laurent expansion in e = (4 — D)/2. 



same colour representation, which is achieved by setting Ca = Cp = 2Tp and rij = 1 in 
the QCD result [27]. It turns out that with this adjustment the leading transcendentality 
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pieces of the quark and gluon form factor become equal at one, two, and three loops in all 
coefficients up to transcendental weight eight, i.e. 0{e^), ©(e^), and 0{e^) at one, two, and 
three loops, respectively. Moreover, the leading transcendentality pieces of the quark and 
gluon form factor coincide - up to a factor of 2^ [L is the number of loops) which is due to 
normalisation - with the coefficients of the scalar form factor in A/" = 4 SYM computed in 
the present work. This again holds true at one, two, and three loops and for all coefficients 
up to weight eight, and serves as an important check of our result. 

The question arises if the leading transcendentality principle [54-56] between QCD 
and A/" = 4 SYM carries over to more general quantites like scattering amplitudes, or if it 
is a special feature of form factors since they have only two external partons. 

In fact, there are counterexamples in the case of scattering amplitudes [19]. For in- 
stance, the = 1 supersymmetric one-loop four-point amplitudes [57] have a leading 
transcendentality piece which is not of the = 4 SYM form, because it has 1 /u power-law 
factors. This makes the property we have found for the form factor even more surprising. 

6. Logarithm of the form factor 

The logarithm of the form factor is given by 



In [Fs) = In 1 + a F'^'' + « ^ 4 + « ^ 4 + ^'(a*) 



+0{a^) , (6.1) 

where 



^ e2 12 3 1440 V 5 36 

4 / 9497r6 _ 49^\ 5 / 127(7 _ 3297r^C3 _ 3l7r^\ 

^ V 120960 18 y V 7 4320 60 ) 

6 /497r2c| 217C3C5 , 18593^8 X 



+^1^ - ^+ 9676800 J ^'-'^ 

f(2) 1 _ , Cs , ^ /39C5 57r\3\ , 2 /^235C| 26237r6\ 
Fs ) - 24i^ + 4^ + ^^^ 72"j+' 1,^ + "60480"; 

3 I^TStt^ _ 4377r^C5 489(7 \ 
^ ' 96 120 2 / 



\ ^^"^ 10 216 1296 ) ^ ' 



^ ^ ^ 3 V ^ y 1620e2 e V 54 3 / 9 25515 



(6.3) 



1077r^C3 1877r\5 2II8IC7 
^ ' 1620 108 144 



-17- 



2 / 1421 1922C3C5 lOSTvr^Cl 9948077rS\ 

"^^ V 4^ 3 ^ 108 17496000 J 

+0{e^) . (6.4) 

The poles of the logarithm of the form factor have the generic structure [58] 



In {Fs) = Y. 



L=l 

,{L) 



+ O(e0), (6.5) 



4(Le)2 2Le 

with the L-loop cusp 7*^^^ and collinear g^^^ anomalous dimensions [59] given by 

00 

7(a) = 0^7^^) = 4a - 4(20^ + 22(40^ + ©(a^) , (6.6) 

L=l 

goia) = f; a^gj^) = -Csa^ + (^4(5 + y C2C3) a' + O(a^) ■ (6.7) 

We observe that the vanishing of the C'(e'^)-term in the logarithm of the two- loop 
form factor [1] appears to be a coincidence, which does not reproduce at three loops. The 
finite part of the = 4 form factor does therefore not exponentiate, as could have been 
conjectured from the two-loop result. 

7. Ultraviolet divergences in higher dimensions 

Scattering amplitudes and form factors in A/" = 4 super Yang-Mills are ultraviolet (UV) 
finite in four dimensions. It is interesting to ask in what dimension, called critical dimension 
Dc, they first develop UV divergences. This question is of theoretical interest in the context 
of the discussion of possible finiteness of A/" = 8 supergravity, see e.g. [38] and references 
therein. More practically, bounds on the critical UV dimension at a given loop order can 
also be a useful cross-check of computations, or constrain the types of loop integrals that 
can appear. 

There is a bound on the critical dimension based on power counting for supergraphs 
and the background field method. The one-loop case is special due to some technical issue 
with ghosts, but there is a bound for L > 1 loops [60,61], 

D,{L)>4+ '^^^~^\ L>1, (7.1) 



such that for D < Dc the theory is UV finite. The bound (7J) depends on the number M of 
supersymmetries that can be realized off-shell. The maximal amount of supersymmetry can 
be realised using an M = 3 harmonic superspace action for A/" = 4 super Yang-Mills [39]. 
Taking thus AA = 3 in ( |7.lD we have 

Dc{L)>A + ^, L>1. (7.2) 



Equation (7T) is a lower bound for Dc, and in some cases it can be too conserva- 



tive. For example, in the case of scattering amplitudes, studying and excluding potential 
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counterterms bounds on the critical dimension can sometimes be improved, see the re- 
views [62,63]. Investigations of UV properties of four-particle scattering amplitudes have 
shown that their ultraviolet behaviour is better than expected [64]. Their critical dimen- 
sion at two and three loops was shown to be 7 and 6, respectively, suggesting the improved 
bound Dc{L) > 4 + 6/L. The one-loop case is exceptional, but for completeness we note 
that Dc{L = 1) = 8 for the four-particle scattering amplitude. 

We can now study the UV properties for > 4 of the form factor that we have 
computed. There is no statement from Eq. (7^) for the one-loop case, but one can easily 
see that Dc{L = 1) = 6. For the two-loop form factor, the bound ( |7.2D is actually saturated 
since the two-loop form factor develops its first ultraviolet divergence at Dc{L = 2) = 6. 
Moreover, it turns out that in D = 6 — 2e dimensions the leading UV-pole is given 
by the leading UV-pole of the two- loop planar ladder diagram Ei , and that E2 has only a 
simple 1/e pole. 

At three loops, Eq. ( |7.2| ) becomes Dc > 16/3. First of all, we see by power counting 
that diagrams F-j and Fio (see Fig. ^ both have a UV divergence in D = 14/3 dimensions, 
which would violate the supersymmetry bound ( |7.2| ). This comes close to explaining why 
their coefficients are zero, and why other integrals having seven or fewer propagators do 
not appear. A small caveat is that it may not always be possible to write the answer in 
a form such that the UV properties are manifest: one could have a linear combination of 
integrals that individually have worse UV properties than expected, but with appropriate 
UV behaviour of the linear combination. However, as we will see presently, we can make 
the UV properties of the three-loop form factor completely manifest. 

At two loops we found that the bound from superspace counting was saturated. We 
can ask whether the same happens at three loops, i.e. do we have Dc{L = 3) = 16/3? It 
turns out that the three-loop form factor is better behaved in the UV than suggested by 
this equation. It is finite in D = 16/3 and only develops a UV divergence at Dc{L = 3) = 6. 
In order to see this, we take the three-loop expression ( 4.13| ) and trade F^, F^ and -F5 for 



the non-UT integrals F| and F^ by means of Eqs. (3^), ( |3.6D , and (3^), respectively, 
which leads to 

p3-ioop ^ ^_^2^) [8 + 2 F3* 2 F4*] - 2 F2 4 F5* - 2 Fg . (7.3) 

Counting numerators as propagators with negative powers, we see that the three integrals 
in the bracket have nine propagators each, whereas the last three integrals have only eight 
propagators. Since there are no sub-divergences in D = 16/3 we can calculate the leading 
UV pole by simply giving all propagators (and also all numerators^) a common mass m 
and by setting the external momenta pi = P2 = 0. Then the first three integrals are finite 
by naive power counting, and the last three integrals become equal, and cancel due to their 
pre-factors. This renders the three- loop form factor finite in D = 16/3 dimensions. One 
can see the UV finiteness of the = 4 SYM form factor in D = 16/3 also in another, more 
elegant way. We start again from Eq. (7.3), and add zero in the disguise of 

+ 2 - 2 F; , (7.4) 

^Whether or not we give a mass to the numerators changes the expressions only by integrals with nine 
propagators each. The latter are finite in D = 16/3 by naiVe power counting. 
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where Fj is F7 (see Fig. ^ with unit numerator. This choice is particularly convenient 
since Fj is a subtopology of both, F2 and F^ . It is obtained from F2 by shrinking the line 
labelled pa in Fig. ^. Alternatively, F^ is obtained from F^ by shrinking line number 7 in 
Fig. ^ In both cases one subsequently has to set the respective numerator to unity. Hence 
we can rewrite (^) as 



fI-^°"p oc (- 



2 (F2 - F;) + 2 (2 F5* - F; - Fg 



[8 Fi + 2 F3* + 2 F4* 
If we adopt for F2 the parametrisation 

{ki , ki+pi, k2, k2+P2, k3-p2, ks+pi, ki+k2, ki 



ks, k2 + k3} 



(7.5) 
(7.6) 



and write (F2 — Fy ) on a common denominator, the numerator of the latter expression 
reads 



kj - {kz-P2y 



(7.7) 



parametrisation ( p.7[ ) for F5 and write (2 F5 — Fy 



and hence vanishes in the aforementioned UV limit. In complete analogy, we take the 

Fg) on a common denominator, whose 

■Vi-V2f \ , (7.8) 



numerator becomes 



{k2 -Pi) 



+ 



{k2-pi) - {k2 

which clearly also vanishes upon taking the UV limit. Hence Eqs. (7.3) and (7.5) make the 
UV properties of the form factor manifest. This is very similar to how the UV properties 
of four-particle amplitudes can be made manifest, see e.g. ref. [64]. 

It is now interesting to investigate the UV properties of the form factor in D = 6 — 2e 
dimensions. Since the vanishing of (F2 — Fj) and (2F5 — Fj — Fg) should be independent 
of the number of dimensions, we can simply look at the expression 

8F1 + 2F3* + 2F4* , (7.9) 

and the corresponding integrals at one and two loops. Introducing a common propagator 
mass and neglecting external momenta one finds 



2D 



UV D=6-2e 
1 — 



^ ^UV D=6^-2e 



8F^ 



UV 



+ 2F3* + 2F4 



=i,UV D=6-2e 



St [ 



21 



-2e 



-3e 



1 


7r2 


e 


f 

6 




1 

+ 2^ 


2^ 



360 

1 7r2 1 
— + — — a$ 



1 

+- 

e 



1 1 

'6^~ 2? 

Cs 13 1 ■ 

— a<T) 

3 12 9 5 



0(6°; 



where Sr is defined in appendix^, and (see e.g. [65]) 



'J>(-i2,i) 



$(-i 2,i) + 



A:=0 



7rln(3) 



[(A; + a)2]«/2 



vrMS) 10 /TTX 



(7.10) 

(7.11) 
(7.12) 
(7.13) 
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and CI2 is the Clausen function. Hence we find that up to three loops the form factor 
at each loop-order has Dc = 6. Moreover, it turns out that for D = 6 — 2e the leading 
UV-pole is at each loop order given by the leading UV-pole of the respective L-loop 
planar ladder diagram. Since at D = 6 — 2e there might be issues due to the presence of 
sub-divergences, we also computed the UV divergences using a different regulator. After 
having taken the soft limit, we re-insert some external momentum into the graph to serve 
as IR regulator, instead of the mass (essentially, one nullifies one of the pi and takes the 
other one off-shell). In this way one obtains massless propagator type integrals which lead 
to the following result 



e2 + 0(e3 



1 5 



1i ^3 

O 



' -I- 2F4 D=6-2e 



-3e 



1 

6^3 



2e 



+ 



1 



4(3 

3 



0{e) 
79 



0(6°) 



(7.14) 



As expected, the leading divergence at L loops is independent of the regulator, while 
the subleading terms are not. However, when considering log(i<5) in the UV limit there 
are only simple 1/e poles up to three loops. Moreover, these poles are identical in both 
regularisation schemes ( [i".10| ) and ( [7.14 ), and read 

^ (4^)3 ■ 
(7.15) 



UVn D=6-2e 



a a 1 
-7+T2 



3 



17 
18 



- (^-^1 +0(a^6°) 



with Q 



Let us now discuss this result. 

Despite the fact that the form factor is better behaved in the UV than expected, 
one may wonder why the four-particle amplitudes at one- and two loops are even better 
behaved in the UV than the form factor. This is due to the fact that there are specific 
counterterms for the local composite operator 0{x) in higher dimensions. Another way of 
saying this is in terms of operator mixing. We note that in D dimensions, the coupling 
constant g has dimension (4 — D)/2. Therefore, in D = G, the operator tr (cp'^) can mix 
at one loop with the operator g'^ Dtr (</>'^), and other operators having the same quantum 
numbers (we have dropped 5'C/(4) indices for simplicity). Another reason for the better 
UV behaviour of the four-point amplitudes, at least in the planar limit, is the fact that 
amplitudes have a dual conformal symmetry, which implies that the difference between 
the number of propagator factors and numerator factors is four for any loop, whereas form 
factors are not dual conformal invariant and therefore can have fewer propagators per loop. 



8. Discussion and conclusion 

In this paper, we extended the calculation of the two-particle form factor in = 4 SYM 
of ref. [1] to the three-loop order. We employed the unitarity-based method to obtain the 
answer in terms of loop integrals. The result contains both planar and non-planar integrals. 
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The form factor can be expressed in several ways in terms of loop integrals that make 
different properties manifest. One way of writing it, Eq. ( |4.13D is in terms of integrals all 
having uniform transcendentality (UT). Other forms, Eqs. and ( [7.5[ ), do not have this 
property, but in turn have the advantage of making the ultraviolet properties of the form 
factor manifest. In order to see the connection between the two representations, we derived 
identities between non-planar integrals based on reparametrisation invariances. 

We evaluated the form factor in dimensional regularisation by reexpressing the integrals 
appearing in it in terms of conventionally used master integrals, c.f. Eq. ( p.ip , whose e 
expansion is known. This allowed us to evaluate the form factor to O(e^). We verified the 
expected exponentiation of infrared divergences, with the correct values at three loops of 
the cusp and collinear anomalous dimensions. 

We observed that the heuristic leading transcendentality principle that relates anoma- 
lous dimensions in QCD with those in = 4 SYM holds also for the form factor. We 
checked this principle to three loops, up to and including terms of transcendental weight 
eight. 

We also studied the ultraviolet (UV) properties of the form factor in higher dimensions. 
We found that at three loops the UV behaviour is better than suggested by a supersym- 
metry argument. Based on power counting one would expect three- loop integrals having 8 
propagators (or nine propagators, and one loop-dependent numerator factor) to diverge in 
D = 16/3 dimensions. However, we find that the particular linear combinations of integrals 
appearing in the form factor is in fact finite in this dimension, and diverges only in D = 6. 
We found a form, Eqs. (|7.3D and (|7.5D , where this is manifest, and computed the leading 
UV divergence of \og{Fs) in = 6 — 2e dimensions. 

There are a number of interesting further directions. 

It is interesting to compare the UV behaviour of the form factor to that of four-particle 
scattering amplitudes. While there are differences due to specific counterterms allowed for 
composite operators, they both share the property of having better UV behaviour than 
expected. It would be interesting if one could understand the UV behaviour of the form 
factor a priori, perhaps based on the absence of potential counterterms, or from string 
theory arguments. 

We remark that the representations of the form factor in terms of UT integrals. 



Eq. ( 4.13 ), or those making its ultraviolet properties manifest, Eq. ( [7.5D , are simpler than 
that in terms of conventionally used master integrals. This may indicate that, even beyond 
M = 4 SYM, there exists a basis of integrals in terms of which the result looks simpler. 
Similar observations about the simplicity of loop integrands and integrals in the case of 
planar scattering amplitudes were also made in refs. [66] and [17]. 

A further extension of this work could be to investigate generalised form factors with 
more on-shell external legs. At one-loop even all-multiplicity results could be envisaged [2- 
5]. At two loops, at least the three-particle form factors should be computable in a relatively 
straightforward manner, since the relevant integrals (two-loop four-point functions with one 
external leg off-shell, [67]) are known from the calculation of QCD amplitudes for the 1—7-3 
decay kinematics [68,69]. 
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The form factor studied in this paper has a very rich structure, similar to that of 
scattering amphtudes. Planar loop integrands of scattering amplitudes, just like tree am- 
plitudes, satisfy powerful recursion relations [66]. It would be extremely interesting to 
extend the applicability of recursion relations to the non-planar case, and the form factor 
studied here is perhaps the simplest case of this type where non-planar integrals appear. 



Acknowledgments 

It is a pleasure to thank L. Dixon for many stimulating discussions and for sharing his in- 
sights on the transcendentality properties of non-planar loop integrals with us. We would 
also like to thank Z. Bern, H. Johansson, S. Naculich, R. Roiban, and E. Sokatchev for 
useful discussions. This project was started during the "Harmony of Scattering Ampli- 
tudes" program at the KITP Santa Barbara, whose hospitality and support we gratefully 
acknowledge. This research was supported in part by the National Science Foundation 
under Grant No. PHY05-51164. The work of TG was supported by the Swiss National 
Science Foundation (SNF) under grant 200020-138206, JMH was supported in part by 
the Department of Energy grant DE-FG02-90ER40542, TH is supported by the Helmholtz 
Alliance "Physics at the Terascale" . Diagrams were drawn with axodraw [70] . 

A. Explicit results of integrals 

In this appendix we list explicit expressions of the integrals that appear as building blocks 
of the form factor. Our integration measure per loop reads 



h 



and we define the pre-factor 

= (4vr)^/2 r(l - e) ■ ^^-'^ 
A generic integral / can be decomposed according to 

/ = [_^2_ .^]n-I,e _^exp^ ^^ 3) 

where L is the number of loops, and the integer n is fixed by dimensional arguments. P'^p 
contains the Laurent expansion about e = 0. 
We start with the one-loop integral 



-)Cxp 



^1 " r(i-2e) • ^^-^^ 



At two loops the integrals read 
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^exp _ r2(l - e)T\e + l)T\-e) 2,T{1 - e)T{2e + l)T\-e) 



^1 



r2(l-2e) 2r(l-3e) 

, 3r(i-2e)r(e + i)r(26 + i)r4(-6) 

+ 4r(l - 3e) ■ ^^-^^ 

An all-order expression for E2 can be found in [71]. The expansion in e reads 

^2 -+?--^-^ ^ + e(8vrC3-117C5)+e ^267C3 - ^ j 

( + 407r2C5 + QCi) + 6^ f-264C5,3 + 2466C3C5 - + 



\^ 10 -,.,0 3Q24 

+0{e^). (A.6) 



At three loops the integrals with uniform transcendentality (UT) are shown in Figs. ^ 
and ^ and read 

pe.p ^ 1 ^ _ ^ _ 3^ , 1 _ 

^ 36e6 12e4 18e3 216e2 e V 6 2 J 

43C| 56 577r6 /2277r^ _ Tvr^ _ ISDCrN 
18~ ~ 68040 ^ ^ V 540 6 3~ J 

The integral F2 is just Ag",^ from [26], 

F,=S|^[-q'-^rJ]-'-''.Fr, 

^2 -+36i^+187 + ^ + 405? + 7V^^ 

137C| 11607r6 / 8297r^C3 TIOtt^Cs 6451(7 \ 
3 5103 V 405 27 9 J 

2/ 1184 lor^n^ ^ 7127r2c| , 5937497^s^^ 



i"— + — 9^ + 1224™ J + ) ■ 

Moreover, we have 

F, = S'r[-<l'-^v]'"~''■Fr, 

3 36e6 ^ 9e4 ^ 9e3 ^ 540.2 e V 3 9 

1352CI ^ 1737r6 ( 2537r''C3 627r2C5 525(7 ^ 



9 1215 V 27 3 2 J 

2^6272 4696C3C5 7127r2c| 1301609^^ , 3 
+e ^^C5,3 ^ ^ ^y^^ ) + 0(0 , (A.9) 
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J 55C3 : 

36e6 12e4 18e3 ^ 

307CI 187977r6 ^ 1497r^C3 , 2397r2C5 21253(7 



]__jr^_55C3 llTT^ 1 /437r\3 599C5 \ 



^ 36e6 12e4 18e3 216e2 ^ e V 6 6 ^ 



18 68040 108 ^ 2 6 

2/8268 5569C3C5 4397r\| 1848737r^ \ 3 
+e ^^C5,3 + g ^QgQOQ- J + 0(6 ) , (A.IO) 

^exp ^ + Z!_ + + + 1 /^ 7l7r\3 ^ 13(5 



^ 12e6 ^ 27e4 ^ 9e3 ^ 540e2 ^ e V 54 

679(1 ^ 399l7r^ ^ / 28377r^C3 ^ 2057r2C5 25135(7 



136080 V 540 9 24 

2/4006^ ^ lOTT^d 141560637rS\ ^^,3, ,,,,, 

+e (53-59(3(5 ^ +0(e^). (A.ll) 

-r ^ 2 S5,3 s3Sb 27 16329600 J 

The integral Fg is just A^^\ from [26], 

i., = 5^[_^2_.^]-2-3..^exp^ 

pexp _ _ l!^!. _ ^ _ 3737r" 1 (Vm^ _ \ 

^6 - +9e6 27e4 9e3 lOSOe^ + , ^ 27 ^^^V 

^169(1 597977r6 ^ ^ /77r^(3 ^ 8507r2(5 18569(7 



9 136080 V 30 9 6 



2/5I88 , 9362(3(5 44367r\| 1078816037r8 \ 3 
(^^(5,3 + ^ 27 81648000 (^"l^) 



Moreover, we have 

F7 = 5F[V-^^]"'"'^•i^7^^ 
pexp ^ 1 ^ _ _ J[l , 1 /'20vr2(3 13(5 



36e6 27e4 9e3 36e2 e V 27 3 

ll 2337r^ /ir;i^4A„ -70^2, 

9 34020 

2 /'248 1244(3(5 1767r2(| ^ 20731l7r8\ 



226(1 2337r^ /iSlvr^ TOtt^ _ 229(7 \ 
+ ^ 1^ 135 + —9— — g— J 



—(53+ """" - "3 ^ +0(e^), (A.13) 

I S5,3 T g 27 20412000 ) ^ > ' ^ > 

^exp ^ + Z!_ _ + ^ + 1 /"^^^S ^27r2(3 



36e6 27e4 9e3 108e2 e V 3 27 
, 98(1 ^ 267r6 ^ ^ /^_!![!^ _ 707r\5 ^ 835(7 



9 8505 V 15 9 6 



2^248 124(3(5 , 5727r2(| 16159^^ \ 3 

+e (^^(5,3 + + - j + 0(0 , (A.14) 
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Figure 8: Diagrams which do not have uniform transcendentality. As before, labels Pa and pb 
on arrow lines indicate an irreducible scalar product {pa + Pb)^ in the respective numerator, and 
diagrams that lack these labels have unit numerator. 
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pexp ^ 1 _ _ ZIC3 _ 437r^ _ 1 /328C5 _ rj^fCA 
9 ^4e6 54e4 9e3 180e2 e V 3 21 ) 

2 295l7r6 / 1021^4^3 GlOvr^Cs , QU%j\ 

- 1701^ - ' V ^35 ^ + ~^J 

/392 _ 11504C3C5 ^ 2876^ _ 85171^X 

V 3 ' 3 27 1020600 / ^ ^ ^ ' 

^exp ^ T{l-efT{-efT{?,e) _ 



12r(l - 4e) 



(A.16) 



The integrals without homogeneous transcendental weight are collected in Fig. p| and read 



T7* oS r 2 -1 — 3— 3e „* exp 

P*exp _ _J_ , ^ 1 /^28C3 ^ 27r2\ ^ /44C3 _ Svr^ ^ TvH 

3 - 9e6 + 27e4 + ^3 I 3 + 9 I + ,2 I 3 9+15 



-26- 



e V 3 27 9 15 y 3 9 9 

^704C3 ^ 1587r6 _ 327r4 _ 1287r2 ^ ^ /_2816C3 2247r2C3 24587r''C3 



567 15 9 V 3 9 135 

1936CI 944C5 _ 2327r\5 _ 2410(7 ^Utt^ 128n^ _ 262tt^\ 
3 "^"3 9 3 9 ~ll> 945~y 



2/ 35152 _ 16082C7 _ 9640C3C5 3527r^C5 3776(5 54167r^C3 5127r 
V 15 ^^'^ 3 3 3 3 27 15" 

7744Cg _ 2247r'^C3 89671^3 11264(3 _ 956008?!^ 10487r^ _ 204871^ ~ 



3 9 9 3 637875 945 9 

+0{e^) , (A.17) 



rt3 r 2 -1 — 3— 3e „* exp 

S'r [-q -Fa 



1 5 1 / 10 TT^X 1 / 55(3 40 77r2 



18e6 18e5 V 9 6 7 e3 V 9 9 9 

( 



^ 1 / 136(3 160 ^ 287r2 llvr^X 1 544(3 4371^(3 599(5 



e2 V 9 9 ' 9 108 y ' e V 9 ' 3 3 

_^640 _ 1127r2 _ 177r'^\ _ 1108(5 _ 307(| ^ 887r2(3 2176(3 187977r6 



9 9 54 / 3 9 9 9 34020 



^347r'' ^ 4487r2 _ 2560 ^ ^ / 8704(3 _ 3527r2(3 _ 1497r^(3 _ 7360(| ^ 4432(5 



27 99V9 9 54 9 3 

^^^2^ 21253(7 10240 17927r2 1367r'i 30557r6\ ^ f IQbm ^ 
+2397r (5 - + ^J+^ {-^^'^^ 

, 11138(3(5 , 2^ 17728(5 4397r\| , 29440(| 48467r\3 
-17273(7 + + 180;r (5 - ^ + — ^ ^ 

14087r^(3 _ 34816(3 _ 1848737r^ 122207r6 5447r^ 71687r^ _ 40960 \ 
^ 9 9 54000 ^ 1701 ^ 27 ^ 9 9~y 

+0{e^) , (A.18) 

_ 1 5 1 /5 7r2\ 1 /57r2 20\ 
"18^6 " 36^5 +^4 1^9 + +^ Vis" ~ yj 

+1 + ^ _ _ + i _ 77^' C3 160(5 

e2 V 9 9 9 40 y e V 9 54 3 

_320 407r^_597r^\ 4^ 

9 9 540 y/ ^^18 9 68040 

135 9 9 V 9 12 9 ^ 

4457r^(5 74815(7 5120 6407r^ _ 2367r^ 25197r6\ 2 ( 12518 
9 24 135~ 2430 J ^ ^ V 5~''^'^ 



6 ' ' 3 '" 27 9 135 

512C3 1 32 4 8 2 577rS 50387r6 9447r^ 25607r2 20480\ 

-1287r2C3 + — + + + 

^ 9 5832000 1215 135 9 9 J 

+0{e^) , (A.19) 



^exp ^ _^ ^ 1 _ 22C3\ ^ 1 /88C3 _ ISvr^ _ ^ II8C5 287r2C3 



«^ 9e3 ^ e2 V 9 3 J e \ 3 9 15 y 3 9 

352C3 _^ IGTr"^ _^ 647r2 _^ ^ / 1408(3 _^ 1127r2C3 _^ 968C| 472(5 256n' 



315 9 V3 9 3 3 9 

647r4 ^ 13l7r^\ ^ ^2 5632(3 _ 4487r2C3 ^ 1127r^C3 _ 3872C| ^ I888C5 



1767r2C5 8041C7 10247r2 2567r^ 5247r6 \ ^^,0, , 
+ ^ + ^7^ + — ^ — + — TT (A.20) 



15 945 7 V 3 9 9 3 3 

567r^ _ 5247r^ 
I5 945" 

^exp ^ _ J_ + A + f ^ + ^ _ + 1 ^ 272C3 320 



«2 36e5 9e4^eHl8 9^ 62^9 9 9 ^ 6^ 9 9 
^ 567r2 ^ ITtt^X ^ 554(5 _ 447r2C3 ^ IO88C3 _ ITtt^ _ 2247r2 ^ 1280 



108 y 3 9 9 27 9 9 

/1767r2C3 _ 4352(3 3680C| 22I6C5 _ 5120 8967r2 GSn* 30557r6\ 
'^^ \ 9 9 9 3 9~ 9 ^ 27 3402 j 

2 frmSCs 7047r2C3 ^ 24237r4C3 14720C| ^ 8864(5 2. ^ 17273(7 

1—9 9— + ^35 9— + + 

20480 35847r2 2727r^ 61107r6\ . 3 



+ 



9 9 27 1701 



-j+0{e^). (A.21) 

At three loops we also cross-checked the major part of the integrals with the sector 
decomposition program FIESTA [72,73]. 



B. Form factor in terms of master integrals 

Just as in QCD, the three-loop scalar form factor in A/" = 4 can be reduced to master inte- 
grals by means of the Laporta algorithm [44] , for which we used the program REDUZE [45] . 
One obtains 



.(3) _ ;,3 r , (3^ - , 
S r(£»-4)(5L>-22)^''' 

20(3D - 13) (D - 3) 



2(3L»-14) 4(2Ll-9)(3i:>- 14) 

Ag o — A; 



{D-A){2D-9) 
16(3D - 13){3D - 11) 
(2D-9)(5L>-22) 



7,1 



5D-22 
40{D - 3) 



D 



A7,2 + 



(D-4)(5£)-22) 
8{D - 4) 



8,1 



(2L>-9)(5L>-22) 



A' 



7,3 



^ 16(3Z^ - 13)(3Z^ - 11) 
^'^ (2L>-9)(5L>-22) ^'^ 
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A, 



6,2 



128(2L> - 7)(L> -3)2 

'3(L» - 4){3D - U){5D - 22) ^'^ 

16(2Ll - 7)(5L> - 18) (521)2 _ 435^ _^ ^^28) 

9(L>-4)2(2L>-9)(5L>-22) 

16(2£' - 7)(3£' - 14)(3£> - 10) - 3) 

(L>- 4)3(51) -22) 

128(2i:> - 7)(3D - 8) (911)2 _ 3211) + I85l) {D - 3)2 

3(D-4)4(2D-9)(5D-22) 

128(21) - 7) (14971)3 - 204231)2 + 92824Z) - 140556) {D - 3f 

^5,2 



A, 



6,3 



5,1 



9(Z) - 4)4(21) - 9)(3Z) - 14)(5D - 22) 

' , 48(3Z)-10)(D- 
-^6,1 + _ 4)3 

16(3D - 10)(3Z) - 8) (144Z)2 - 1285D + 2866) (D - 3)2 



, 4(D - 3) ^ 64(1) -3)3 ^ , 48(3Z)-10)(D-3)2 ^ 
+"D^ + (D-4)3 + CD^rip ^^'^ 



+ 



(D - 4)4(21) - 9)(5D - 22) 
1771)2 - 1584D + 3! 
-4)4(2D -9)(5D - 
64(2Z)-5)(3D-8)(Z)-3) 



^5, 



128(2i:> - 7) (1771)2 - 1584D + 3542) (D - 3)3 
^ 3(D-4)4(2D-9)(5D-22) ^^'^ 



9(D - 4)5(2Z) - 9)(3Z) - 14)(5Z) - 22) 

X (25021)^ - 51273D^ + 419539D^ - 17136881)2 + 3495112Z) - 2848104) ^4,1 



2 



4(i)-3) , 48(3Z)-10)(D-3) ^ 



(B.l) 



-Re is given in Eq. ( p.l4 ). In order to arrive at Eq. ( |5.2| ) we have to plug in D = 4 — 2e and 



the e-expansions for the master integrals from Eqs. (A. 7) - (A. 27) of [30], together with 
their higher order e-terms from [32]. 

C. Four-point amplitude to two loops 

Here we summarise the known four-point amplitude in A/" = 4 super Yang-Mills to two 
loop order. As we have seen in the main text, both leading and subleading terms in colour 
are required when computing the form factor at leading colour using unitarity. 

We consider four-point amplitudes in SU{N) gauge theories with all particles in the 
adjoint representation. Let us review the decomposition of the latter into a trace basis 
with partial amplitudes as coefficients [74,75]. 

At tree-level, we have 

At""" = g"^ fi^' ^ TiiT''-WT''-(^)T''-(3)T''-W)A%l^^{ail),ai2),ai3),ai4:)), (C.l) 
(Te54/Z4 

where sum goes over the six non-cyclic permutations of (1234), i.e. S^/Z^ = {(1234), (2134), 
(1243), (2314), (3241), (3214)}. The AX^^r^ are 'partial amplitudes'. The arguments of ^ 
and A in Eq. ( p.l]) are abbreviations, i.e. 1 stands for a given particle (gluon, fermion, or 
scalar) of a given helicity and momentum p'^. The T'* are the (A'^2 _ 1) matrices in the 
fundamental representation of SU{N). 



-29- 




Figure 9: Scalar box integrals appearing in four-particle amplitudes to two loops 

At loop level, double trace terms are present as well. Other possible trace terms vanish 
since Tr(T") = for SU{N). We have, at one loop 

ji'ioop ^ gA^ie ArTr(r'^-(i)r'^-(2)r'^-{3)r''-W)yl^;-^;°°f((7(l),(7(2),a(3),a(4)) 

erg 54 /•^4 

+5V^' Yl Tr(T'^-(i)T'^-(2))Tr(T'^-(3)rM4))Ai-^°°P(a(l),f7(2),CT(3),(T(4)), (C.2) 
and two loops [53], 

^X^-ioop ^ g6^6e Tr(T'*-(i)T'*-(2)r''-{3)r''-w) X (C.3) 

crGS4/Z4 

x[N'Al-[r''^''{a{l),a{2),a{3),a{A))+Al-[r^^^ 

+9^1^^' A^Tr(r"-{i)r''-{2))Tr(r''-{3)r'*-(4))^2-/oop^^^-^^^^^2),cj(3),(T(4)) . 

Here S^/Z^ = {(1234), (1324), (1423)}. The double trace terms are subleading in the 
expansion in powers of N. At the two-loop order, we also have the appearence of subleading- 
in-A^ terms in the single trace terms, denoted by the superscript SC, while the leading-in-A^ 
terms have superscript LC. 

N = 4: super symmetric Ward identities imply that for MHV amplitudes the loop-level 
amplitudes are proportional to the tree-level ones, for any choice of external particles and 
helicities. We have 

Al-;r(l,2,3,4) = -st4;ffi(l,2,3,4)xl-'"°P(s,t), (C.4) 

where^ 



fc2(A: - piY{k -pi- p2Y{k + p^Y ' ^^-^^ 

is the one-loop scalar box integral, see Fig. H The remaining subleading colour amplitudes 
at one loop are all equal and given by 

A\;i7= E vll-i;r(a(l),a(2),a(3),a(4)), (C.6) 



Note that our convention of defining loop integrals differs from that of ref. [53] by a factor of i per loop 



order, cf Eq. (A.l) 
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which is the consequence of a U{1) decouphng identity [74]. 

At two loops, the partial amplitudes leading in are given by [53] 

^Sr'^^(l>2,3,4) = +st4':i^fi(l,2,3,4) {sl[{s,t) + stl[{t,s)) , (C.7) 

where l[{s,t) is the planar double box integral, see Fig. ^ 
The partial amplitudes subleading in N are given by [53] 

j^2-^ioop,sc ^-^^ 2, 3, 4) = 2ylf (1, 2; 3, 4) + 2Af (3, 4; 2, 1) + 2A^{1, 4; 2, 3) + 2Af (2, 3; 4, 1) 

- 4<(1, 3; 2, 4) - 4^(2, 4; 3, 1) + 2Af ^(1; 2; 3, 4) + 2^1^ ^(3; 4; 2, 1) 

+ 2<^(1; 4; 2, 3) + 2<^(2; 3; 4, 1) - 4<^(1; 3; 2, 4) - 4<^(2; 4; 3, 1) , (C.8) 

and 

Al-fP{l; 2; 3, 4) = 6ylf (1, 2; 3, 4) + 6Af (1, 2; 4, 3) + 4<^(1; 2; 3, 4) + 4<^(3; 4; 2, 1) 

- 2<^(1; 4; 2, 3) - 2<^(2; 3; 4, 1) - 2<^(1; 3; 2, 4) - 2<^(2; 4; 3, 1) , (C.9) 

where 

^f(l,2;3,4) =s22S23 4':ri(l' 2,3,4) Xf(si2,S23), (C.IO) 

2; 3, 4) ^ s^^^ss ^-iKl, 2, 3, 4)Xf ^(.12, S23) , (C.ll) 

and where and X^-^ are the planar and non-planar double box integral, respectively, 
see Fig. |[ 

We remark that the expression for the double trace terms A^7i'^°^ can be obtained 
from the single trace terms using identities derived from group theory [75,76]. 
The tree-level amplitude we need has external scalars only. It is given by 

4':iri('^12(l)><^12(2),<^34(3),<A34(4)) = -i^. (C.12) 
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